Incorporating heterogeneity in farmer disease control behaviour into a livestock disease transmission model.

Human behaviour is critical to effective responses to livestock disease outbreaks, especially with respect to vaccination uptake. Traditionally, mathematical models used to inform this behaviour have not taken heterogeneity in farmer behaviour into account. We address this by exploring how heterogeneity in farmers vaccination behaviour can be incorporated to inform mathematical models. We developed and used a graphical user interface to elicit farmers (n = 60) vaccination decisions to an unfolding fast-spreading epidemic and linked this to their psychosocial and behavioural profiles. We identified, via cluster analysis, robust patterns of heterogeneity in vaccination behaviour. By incorporating these vaccination behavioural groupings into a mathematical model for a fast-spreading livestock infection, using computational simulation we explored how the inclusion of heterogeneity in farmer disease control behaviour may impact epidemiological and economic focused outcomes. When assuming homogeneity in farmer behaviour versus configurations informed by the psychosocial profile cluster estimates, the modelled scenarios revealed a disconnect in projected distributions and threshold statistics across outbreak size, outbreak duration and economic metrics.

Targeting of multiple tumor-associated antigens by individual T cell receptors during successful cancer immunotherapy.

The T cells of the immune system can target tumors and clear solid cancers following tumor-infiltrating lymphocyte (TIL) therapy. We used combinatorial peptide libraries and a proteomic database to reveal the antigen specificities of persistent cancer-specific T cell receptors (TCRs) following successful TIL therapy for stage IV malignant melanoma. Remarkably, individual TCRs could target multiple different tumor types via the HLA A∗02:01-restricted epitopes EAAGIGILTV, LLLGIGILVL, and NLSALGIFST from Melan A, BST2, and IMP2, respectively. Atomic structures of a TCR bound to all three antigens revealed the importance of the shared x-x-x-A/G-I/L-G-I-x-x-x recognition motif. Multi-epitope targeting allows individual T cells to attack cancer in several ways simultaneously. Such "multipronged" T cells exhibited superior recognition of cancer cells compared with conventional T cell recognition of individual epitopes, making them attractive candidates for the development of future immunotherapies.

Health economic evaluation of strategies to eliminate gambiense human African trypanosomiasis in the Mandoul disease focus of Chad.

Human African trypanosomiasis, caused by the gambiense subspecies of Trypanosoma brucei (gHAT), is a deadly parasitic disease transmitted by tsetse. Partners worldwide have stepped up efforts to eliminate the disease, and the Chadian government has focused on the previously high-prevalence setting of Mandoul. In this study, we evaluate the economic efficiency of the intensified strategy that was put in place in 2014 aimed at interrupting the transmission of gHAT, and we make recommendations on the best way forward based on both epidemiological projections and cost-effectiveness. In our analysis, we use a dynamic transmission model fit to epidemiological data from Mandoul to evaluate the cost-effectiveness of combinations of active screening, improved passive screening (defined as an expansion of the number of health posts capable of screening for gHAT), and vector control activities (the deployment of Tiny Targets to control the tsetse vector). For cost-effectiveness analyses, our primary outcome is disease burden, denominated in disability-adjusted life-years (DALYs), and costs, denominated in 2020 US$. Although active and passive screening have enabled more rapid diagnosis and accessible treatment in Mandoul, the addition of vector control provided good value-for-money (at less than $750/DALY averted) which substantially increased the probability of reaching the 2030 elimination target for gHAT as set by the World Health Organization. Our transmission modelling and economic evaluation suggest that the gains that have been made could be maintained by passive screening. Our analysis speaks to comparative efficiency, and it does not take into account all possible considerations; for instance, any cessation of ongoing active screening should first consider that substantial surveillance activities will be critical to verify the elimination of transmission and to protect against the possible importation of infection from neighbouring endemic foci.

Modelling to infer the role of animals in gambiense human African trypanosomiasis transmission and elimination in the DRC.

Gambiense human African trypanosomiasis (gHAT) has been targeted for elimination of transmission (EoT) to humans by 2030. Whilst this ambitious goal is rapidly approaching, there remain fundamental questions about the presence of non-human animal transmission cycles and their potential role in slowing progress towards, or even preventing, EoT. In this study we focus on the country with the most gHAT disease burden, the Democratic Republic of Congo (DRC), and use mathematical modelling to assess whether animals may contribute to transmission in specific regions, and if so, how their presence could impact the likelihood and timing of EoT. By fitting two model variants-one with, and one without animal transmission-to the human case data from 2000-2016 we estimate model parameters for 158 endemic health zones of the DRC. We evaluate the statistical support for each model variant in each health zone and infer the contribution of animals to overall transmission and how this could impact predicted time to EoT. We conclude that there are 24/158 health zones where there is substantial to decisive statistical support for some animal transmission. However-even in these regions-we estimate that animals would be extremely unlikely to maintain transmission on their own. Animal transmission could hamper progress towards EoT in some settings, with projections under continuing interventions indicating that the number of health zones expected to achieve EoT by 2030 reduces from 68/158 to 61/158 if animal transmission is included in the model. With supplementary vector control (at a modest 60% tsetse reduction) added to medical screening and treatment interventions, the predicted number of health zones meeting the goal increases to 147/158 for the model including animal transmission. This is due to the impact of vector reduction on transmission to and from all hosts.

Identifying regions for enhanced control of gambiense sleeping sickness in the Democratic Republic of Congo.

Gambiense human African trypanosomiasis (sleeping sickness, gHAT) is a disease targeted for elimination of transmission by 2030. While annual new cases are at a historical minimum, the likelihood of achieving the target is unknown. We utilised modelling to study the impacts of four strategies using currently available interventions, including active and passive screening and vector control, on disease burden and transmission across 168 endemic health zones in the Democratic Republic of the Congo. Median projected years of elimination of transmission show only 98 health zones are on track despite significant reduction in disease burden under medical-only strategies (64 health zones if > 90% certainty required). Blanket coverage with vector control is impractical, but is predicted to reach the target in all heath zones. Utilising projected disease burden under the uniform medical-only strategy, we provide a priority list of health zones for consideration for supplementary vector control alongside medical interventions.

Cost-effectiveness of sleeping sickness elimination campaigns in five settings of the Democratic Republic of Congo.

Gambiense human African trypanosomiasis (gHAT) is marked for elimination of transmission by 2030, but the disease persists in several low-income countries. We couple transmission and health outcomes models to examine the cost-effectiveness of four gHAT elimination strategies in five settings - spanning low- to high-risk - of the Democratic Republic of Congo. Alongside passive screening in fixed health facilities, the strategies include active screening at average or intensified coverage levels, alone or with vector control with a scale-back algorithm when no cases are reported for three consecutive years. In high or moderate-risk settings, costs of gHAT strategies are primarily driven by active screening and, if used, vector control. Due to the cessation of active screening and vector control, most investments (75-80%) are made by 2030 and vector control might be cost-saving while ensuring elimination of transmission. In low-risk settings, costs are driven by passive screening, and minimum-cost strategies consisting of active screening and passive screening lead to elimination of transmission by 2030 with high probability.

Update of transmission modelling and projections of gambiense human African trypanosomiasis in the Mandoul focus, Chad.

BACKGROUND: In recent years, a programme of vector control, screening and treatment of gambiense human African trypanosomiasis (gHAT) infections led to a rapid decline in cases in the Mandoul focus of Chad. To represent the biology of transmission between humans and tsetse, we previously developed a mechanistic transmission model, fitted to data between 2000 and 2013 which suggested that transmission was interrupted by 2015. The present study outlines refinements to the model to: (1) Assess whether elimination of transmission has already been achieved despite low-level case reporting; (2) quantify the role of intensified interventions in transmission reduction; and (3) predict the trajectory of gHAT in Mandoul for the next decade under different strategies. METHOD: Our previous gHAT transmission model for Mandoul was updated using human case data (2000-2019) and a series of model refinements. These include how diagnostic specificity is incorporated into the model and improvements to the fitting method (increased variance in observed case reporting and how underreporting and improvements to passive screening are captured). A side-by-side comparison of fitting to case data was performed between the models. RESULTS: We estimated that passive detection rates have increased due to improvements in diagnostic availability in fixed health facilities since 2015, by 2.1-fold for stage 1 detection, and 1.5-fold for stage 2. We find that whilst the diagnostic algorithm for active screening is estimated to be highly specific (95% credible interval (CI) 99.9-100%, Specificity = 99.9%), the high screening and low infection levels mean that some recently reported cases with no parasitological confirmation might be false positives. We also find that the focus-wide tsetse reduction estimated through model fitting (95% CI 96.1-99.6%, Reduction = 99.1%) is comparable to the reduction previously measured by the decline in tsetse catches from monitoring traps. In line with previous results, the model suggests that transmission was interrupted in 2015 due to intensified interventions. CONCLUSIONS: We recommend that additional confirmatory testing is performed in Mandoul to ensure the endgame can be carefully monitored. More specific measurement of cases, would better inform when it is safe to stop active screening and vector control, provided there is a strong passive surveillance system in place.

Synthetic Peptides with Inadvertent Chemical Modifications Can Activate Potentially Autoreactive T Cells.

The human CD8+ T cell clone 6C5 has previously been shown to recognize the tert-butyl-modified Bax161-170 peptide LLSY(3-tBu)FGTPT presented by HLA-A*02:01. This nonnatural epitope was likely created as a by-product of fluorenylmethoxycarbonyl protecting group peptide synthesis and bound poorly to HLA-A*02:01. In this study, we used a systematic approach to identify and characterize natural ligands for the 6C5 TCR. Functional analyses revealed that 6C5 T cells only recognized the LLSYFGTPT peptide when tBu was added to the tyrosine residue and did not recognize the LLSYFGTPT peptide modified with larger (di-tBu) or smaller chemical groups (Me). Combinatorial peptide library screening further showed that 6C5 T cells recognized a series of self-derived peptides with dissimilar amino acid sequences to LLSY(3-tBu)FGTPT. Structural studies of LLSY(3-tBu)FGTPT and two other activating nonamers (IIGWMWIPV and LLGWVFAQV) in complex with HLA-A*02:01 demonstrated similar overall peptide conformations and highlighted the importance of the position (P) 4 residue for T cell recognition, particularly the capacity of the bulky amino acid tryptophan to substitute for the tBu-modified tyrosine residue in conjunction with other changes at P5 and P6. Collectively, these results indicated that chemical modifications directly altered the immunogenicity of a synthetic peptide via molecular mimicry, leading to the inadvertent activation of a T cell clone with unexpected and potentially autoreactive specificities.

Viral Molecular Mimicry Influences the Antitumor Immune Response in Murine and Human Melanoma.

Molecular mimicry is one of the leading mechanisms by which infectious agents can induce autoimmunity. Whether a similar mechanism triggers an antitumor immune response is unexplored, and the role of antiviral T cells infiltrating the tumor has remained anecdotal. To address these questions, we first developed a bioinformatic tool to identify tumor peptides with high similarity to viral epitopes. Using peptides identified by this tool, we demonstrated that, in mice, preexisting immunity toward specific viral epitopes enhanced the efficacy of cancer immunotherapy via molecular mimicry in different settings. To understand whether this mechanism could partly explain immunotherapy responsiveness in humans, we analyzed a cohort of patients with melanoma undergoing anti-PD1 treatment who had a high IgG titer for cytomegalovirus (CMV). In this cohort of patients, we showed that high levels of CMV-specific antibodies were associated with prolonged progression-free survival and found that, in some cases, peripheral blood mononuclear cells (PBMC) could cross-react with both melanoma and CMV homologous peptides. Finally, T-cell receptor sequencing revealed expansion of the same CD8+ T-cell clones when PBMCs were expanded with tumor or homologous viral peptides. In conclusion, we have demonstrated that preexisting immunity and molecular mimicry could influence the response to immunotherapies. In addition, we have developed a free online tool that can identify tumor antigens and neoantigens highly similar to pathogen antigens to exploit molecular mimicry and cross-reactive T cells in cancer vaccine development.

Quantifying epidemiological drivers of gambiense human African Trypanosomiasis across the Democratic Republic of Congo.

Gambiense human African trypanosomiasis (gHAT) is a virulent disease declining in burden but still endemic in West and Central Africa. Although it is targeted for elimination of transmission by 2030, there remain numerous questions about the drivers of infection and how these vary geographically. In this study we focus on the Democratic Republic of Congo (DRC), which accounted for 84% of the global case burden in 2016, to explore changes in transmission across the country and elucidate factors which may have contributed to the persistence of disease or success of interventions in different regions. We present a Bayesian fitting methodology, applied to 168 endemic health zones (∼100,000 population size), which allows for calibration of a mechanistic gHAT model to case data (from the World Health Organization HAT Atlas) in an adaptive and automated framework. It was found that the model needed to capture improvements in passive detection to match observed trends in the data within former Bandundu and Bas Congo provinces indicating these regions have substantially reduced time to detection. Health zones in these provinces generally had longer burn-in periods during fitting due to additional model parameters. Posterior probability distributions were found for a range of fitted parameters in each health zone; these included the basic reproduction number estimates for pre-1998 (R0) which was inferred to be between 1 and 1.14, in line with previous gHAT estimates, with higher median values typically in health zones with more case reporting in the 2000s. Previously, it was not clear whether a fall in active case finding in the period contributed to the declining case numbers. The modelling here accounts for variable screening and suggests that underlying transmission has also reduced greatly-on average 96% in former Equateur, 93% in former Bas Congo and 89% in former Bandundu-Equateur and Bandundu having had the highest case burdens in 2000. This analysis also sets out a framework to enable future predictions for the country.

Regulation of Cell Type-Specific Immunity Networks in Arabidopsis Roots.

While root diseases are among the most devastating stresses in global crop production, our understanding of root immunity is still limited relative to our knowledge of immune responses in leaves. Considering that root performance is based on the concerted functions of its different cell types, we undertook a cell type-specific transcriptome analysis to identify gene networks activated in epidermis, cortex, and pericycle cells of Arabidopsis (Arabidopsis thaliana) roots challenged with two immunity elicitors, the bacterial flagellin-derived flg22 and the endogenous Pep1 peptide. Our analyses revealed distinct immunity gene networks in each cell type. To further substantiate our understanding of regulatory patterns underlying these cell type-specific immunity networks, we developed a tool to analyze paired transcription factor binding motifs in the promoters of cell type-specific genes. Our study points toward a connection between cell identity and cell type-specific immunity networks that might guide cell types in launching immune response according to the functional capabilities of each cell type.

GPU-Accelerated Discovery of Pathogen-Derived Molecular Mimics of a T-Cell Insulin Epitope.

The strong links between (Human Leukocyte Antigen) HLA, infection and autoimmunity combine to implicate T-cells as primary triggers of autoimmune disease (AD). T-cell crossreactivity between microbially-derived peptides and self-peptides has been shown to break tolerance and trigger AD in experimental animal models. Detailed examination of the potential for T-cell crossreactivity to trigger human AD will require means of predicting which peptides might be recognised by autoimmune T-cell receptors (TCRs). Recent developments in high throughput sequencing and bioinformatics mean that it is now possible to link individual TCRs to specific pathologies for the first time. Deconvolution of TCR function requires knowledge of TCR specificity. Positional Scanning Combinatorial Peptide Libraries (PS-CPLs) can be used to predict HLA-restriction and define antigenic peptides derived from self and pathogen proteins. In silico search of the known terrestrial proteome with a prediction algorithm that ranks potential antigens in order of recognition likelihood requires complex, large-scale computations over several days that are infeasible on a personal computer. We decreased the time required for peptide searching to under 30 min using multiple blocks on graphics processing units (GPUs). This time-efficient, cost-effective hardware accelerator was used to screen bacterial and fungal human pathogens for peptide sequences predicted to activate a T-cell clone, InsB4, that was isolated from a patient with type 1 diabetes and recognised the insulin B-derived epitope HLVEALYLV in the context of disease-risk allele HLA A*0201. InsB4 was shown to kill HLA A*0201+ human insulin producing ß-cells demonstrating that T-cells with this specificity might contribute to disease. The GPU-accelerated algorithm and multispecies pathogen proteomic databases were validated to discover pathogen-derived peptide sequences that acted as super-agonists for the InsB4 T-cell clone. Peptide-MHC tetramer binding and surface plasmon resonance were used to confirm that the InsB4 TCR bound to the highest-ranked peptide agonists derived from infectious bacteria and fungi. Adoption of GPU-accelerated prediction of T-cell agonists has the capacity to revolutionise our understanding of AD by identifying potential targets for autoimmune T-cells. This approach has further potential for dissecting T-cell responses to infectious disease and cancer.

PeTTSy: a computational tool for perturbation analysis of complex systems biology models.

BACKGROUND: Over the last decade sensitivity analysis techniques have been shown to be very useful to analyse complex and high dimensional Systems Biology models. However, many of the currently available toolboxes have either used parameter sampling, been focused on a restricted set of model observables of interest, studied optimisation of a objective function, or have not dealt with multiple simultaneous model parameter changes where the changes can be permanent or temporary. RESULTS: Here we introduce our new, freely downloadable toolbox, PeTTSy (Perturbation Theory Toolbox for Systems). PeTTSy is a package for MATLAB which implements a wide array of techniques for the perturbation theory and sensitivity analysis of large and complex ordinary differential equation (ODE) based models. PeTTSy is a comprehensive modelling framework that introduces a number of new approaches and that fully addresses analysis of oscillatory systems. It examines sensitivity analysis of the models to perturbations of parameters, where the perturbation timing, strength, length and overall shape can be controlled by the user. This can be done in a system-global setting, namely, the user can determine how many parameters to perturb, by how much and for how long. PeTTSy also offers the user the ability to explore the effect of the parameter perturbations on many different types of outputs: period, phase (timing of peak) and model solutions. PeTTSy can be employed on a wide range of mathematical models including free-running and forced oscillators and signalling systems. To enable experimental optimisation using the Fisher Information Matrix it efficiently allows one to combine multiple variants of a model (i.e. a model with multiple experimental conditions) in order to determine the value of new experiments. It is especially useful in the analysis of large and complex models involving many variables and parameters. CONCLUSIONS: PeTTSy is a comprehensive tool for analysing large and complex models of regulatory and signalling systems. It allows for simulation and analysis of models under a variety of environmental conditions and for experimental optimisation of complex combined experiments. With its unique set of tools it makes a valuable addition to the current library of sensitivity analysis toolboxes. We believe that this software will be of great use to the wider biological, systems biology and modelling communities.

Identification of human viral protein-derived ligands recognized by individual MHCI-restricted T-cell receptors.

Evidence indicates that autoimmunity can be triggered by virus-specific CD8(+) T cells that crossreact with self-derived peptide epitopes presented on the cell surface by major histocompatibility complex class I (MHCI) molecules. Identification of the associated viral pathogens is challenging because individual T-cell receptors can potentially recognize up to a million different peptides. Here, we generate peptide length-matched combinatorial peptide library (CPL) scan data for a panel of virus-specific CD8(+) T-cell clones spanning different restriction elements and a range of epitope lengths. CPL scan data drove a protein database search limited to viruses that infect humans. Peptide sequences were ranked in order of likelihood of recognition. For all anti-viral CD8(+) T-cell clones examined in this study, the index peptide was either the top-ranked sequence or ranked as one of the most likely sequences to be recognized. Thus, we demonstrate that anti-viral CD8(+) T-cell clones are highly focused on their index peptide sequence and that 'CPL-driven database searching' can be used to identify the inciting virus-derived epitope for a given CD8(+) T-cell clone. Moreover, to augment access to CPL-driven database searching, we have created a publicly accessible webtool. Application of these methodologies in the clinical setting may clarify the role of viral pathogens in the etiology of autoimmune diseases.

CSI: a nonparametric Bayesian approach to network inference from multiple perturbed time series gene expression data.

Here we introduce the causal structure identification (CSI) package, a Gaussian process based approach to inferring gene regulatory networks (GRNs) from multiple time series data. The standard CSI approach infers a single GRN via joint learning from multiple time series datasets; the hierarchical approach (HCSI) infers a separate GRN for each dataset, albeit with the networks constrained to favor similar structures, allowing for the identification of context specific networks. The software is implemented in MATLAB and includes a graphical user interface (GUI) for user friendly inference. Finally the GUI can be connected to high performance computer clusters to facilitate analysis of large genomic datasets.

Modeling meiotic chromosomes indicates a size dependent contribution of telomere clustering and chromosome rigidity to homologue juxtaposition.

Meiosis is the cell division that halves the genetic component of diploid cells to form gametes or spores. To achieve this, meiotic cells undergo a radical spatial reorganisation of chromosomes. This reorganisation is a prerequisite for the pairing of parental homologous chromosomes and the reductional division, which halves the number of chromosomes in daughter cells. Of particular note is the change from a centromere clustered layout (Rabl configuration) to a telomere clustered conformation (bouquet stage). The contribution of the bouquet structure to homologous chromosome pairing is uncertain. We have developed a new in silico model to represent the chromosomes of Saccharomyces cerevisiae in space, based on a worm-like chain model constrained by attachment to the nuclear envelope and clustering forces. We have asked how these constraints could influence chromosome layout, with particular regard to the juxtaposition of homologous chromosomes and potential nonallelic, ectopic, interactions. The data support the view that the bouquet may be sufficient to bring short chromosomes together, but the contribution to long chromosomes is less. We also find that persistence length is critical to how much influence the bouquet structure could have, both on pairing of homologues and avoiding contacts with heterologues. This work represents an important development in computer modeling of chromosomes, and suggests new explanations for why elucidating the functional significance of the bouquet by genetics has been so difficult.

Quantitative analysis of regulatory flexibility under changing environmental conditions.

The circadian clock controls 24-h rhythms in many biological processes, allowing appropriate timing of biological rhythms relative to dawn and dusk. Known clock circuits include multiple, interlocked feedback loops. Theory suggested that multiple loops contribute the flexibility for molecular rhythms to track multiple phases of the external cycle. Clear dawn- and dusk-tracking rhythms illustrate the flexibility of timing in Ipomoea nil. Molecular clock components in Arabidopsis thaliana showed complex, photoperiod-dependent regulation, which was analysed by comparison with three contrasting models. A simple, quantitative measure, Dusk Sensitivity, was introduced to compare the behaviour of clock models with varying loop complexity. Evening-expressed clock genes showed photoperiod-dependent dusk sensitivity, as predicted by the three-loop model, whereas the one- and two-loop models tracked dawn and dusk, respectively. Output genes for starch degradation achieved dusk-tracking expression through light regulation, rather than a dusk-tracking rhythm. Model analysis predicted which biochemical processes could be manipulated to extend dusk tracking. Our results reveal how an operating principle of biological regulators applies specifically to the plant circadian clock.

Robustness from flexibility in the fungal circadian clock.

BACKGROUND: Robustness is a central property of living systems, enabling function to be maintained against environmental perturbations. A key challenge is to identify the structures in biological circuits that confer system-level properties such as robustness. Circadian clocks allow organisms to adapt to the predictable changes of the 24-hour day/night cycle by generating endogenous rhythms that can be entrained to the external cycle. In all organisms, the clock circuits typically comprise multiple interlocked feedback loops controlling the rhythmic expression of key genes. Previously, we showed that such architectures increase the flexibility of the clock's rhythmic behaviour. We now test the relationship between flexibility and robustness, using a mathematical model of the circuit controlling conidiation in the fungus Neurospora crassa. RESULTS: The circuit modelled in this work consists of a central negative feedback loop, in which the frequency (frq) gene inhibits its transcriptional activator white collar-1 (wc-1), interlocked with a positive feedback loop in which FRQ protein upregulates WC-1 production. Importantly, our model reproduces the observed entrainment of this circuit under light/dark cycles with varying photoperiod and cycle duration. Our simulations show that whilst the level of frq mRNA is driven directly by the light input, the falling phase of FRQ protein, a molecular correlate of conidiation, maintains a constant phase that is uncoupled from the times of dawn and dusk. The model predicts the behaviour of mutants that uncouple WC-1 production from FRQ's positive feedback, and shows that the positive loop enhances the buffering of conidiation phase against seasonal photoperiod changes. This property is quantified using Kitano's measure for the overall robustness of a regulated system output. Further analysis demonstrates that this functional robustness is a consequence of the greater evolutionary flexibility conferred on the circuit by the interlocking loop structure. CONCLUSIONS: Our model shows that the behaviour of the fungal clock in light-dark cycles can be accounted for by a transcription-translation feedback model of the central FRQ-WC oscillator. More generally, we provide an example of a biological circuit in which greater flexibility yields improved robustness, while also introducing novel sensitivity analysis techniques applicable to a broader range of cellular oscillators.

Prediction of photoperiodic regulators from quantitative gene circuit models.

Photoperiod sensors allow physiological adaptation to the changing seasons. The prevalent hypothesis is that day length perception is mediated through coupling of an endogenous rhythm with an external light signal. Sufficient molecular data are available to test this quantitatively in plants, though not yet in mammals. In Arabidopsis, the clock-regulated genes CONSTANS (CO) and FLAVIN, KELCH, F-BOX (FKF1) and their light-sensitive proteins are thought to form an external coincidence sensor. Here, we model the integration of light and timing information by CO, its target gene FLOWERING LOCUS T (FT), and the circadian clock. Among other predictions, our models show that FKF1 activates FT. We demonstrate experimentally that this effect is independent of the known activation of CO by FKF1, thus we locate a major, novel controller of photoperiodism. External coincidence is part of a complex photoperiod sensor: modeling makes this complexity explicit and may thus contribute to crop improvement.

Luciferases as reporter genes.

Luciferase is the ideal reporter gene to provide temporal and spatial information on promoter activity in Arabidopsis and other eukaryotes; the noninvasive detection of luminescence and short half-life of luciferase activity allow repeated measurements of individual seedlings over several days to assay dynamic changes in gene expression. Transgenic or transiently transformed plants with a luciferase gene under the control of a promoter of interest are required. Detection of the low level of luminescence produced by the luciferase gene in Arabidopsis requires the use of low light detecting charge-coupled device (CCD) cameras or scintillation counters. This chapter contains protocols on assaying and imaging luciferase in vivo and the automation for high-resolution timecourses.